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In  this  paper  a two-parameter  Welbull  model  is  assumed  for  time  to 


failure  (or  number  of  cycles  to  failure) . The  Inverse  power  law  is 


assumed  to  relate  level  of  stress  and  the  Weibull  scale  parameter  with 


the  Weibull  shape  parameter  independent  of  stress  level 


It  is  supposed  that  a life  test  results  in  Weibull  failure-time 


data  applying  to  the  first  r failures  in  a sample  of  size  n.  The  procedure 


derived  herein  allows  one  to  use  these  data  to  obtain  an  approximate 


lower  confidence  bound  on  the  time  of  the  ith  failure  in  a lot  of  size  m 


selected  from  the  same  population  as  the  sample.  The  precision  of  the 


approximation  is  investigated  for  certain  cases  in  which  i = 1,  that  is 


the  first  failure  in  a lot  is  of  interest.  A method  is  given  for  deter- 
IgljPI  mining  whether  the  approximation  is  sufficiently  precise  for  use,  and  an 
example  of  use  of  the  approximation  is  provided. 
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1.  Introduction 


The  following  applies  In  certain  cases  In  which  life  tests  are  per- 
formed on  a sample  at  two  or  more  levels  of  stress  and  Inferences  about  a 
large,  but  finite.  Independent  sample  (lot)  are  to  be  made  for  a nominal 
stress  level  at  which  no  tests  are  performed.  A (^elbull  model  Is  assumed 
for  time  to  failure  and  the  Inverse  power  law  Is  assumed  to  relate  level 
of  stress  and  the  Welbull  scale  parameter  5. 

Actually,  the  approach  given  here  applies  also  when  failure  times  have 
a lognormal  distribution  and  the  data  are  censored.  It  is  applicable  to 
both  Welbull  and  lognormal  data  whenever  some  known  polynomial  function  of 
degree  k relates  a specified  function  of  level  of  stress  and  the  logarithm 
of  time  to  failure  and  when  the  scale  parameter  C (or  a)  of  the  distribution 
of  the  logarithm  of  time  to  failure  (inverse  of  the  Welbull  or  lognormal 
shape  parameter)  Is  Independent  of  stress  level.  See,  for  example,  Charles 
[2]  and  Saunders  [18]. 

The  life  tests  from  which  the  failure-time  data  are  obtained  may  be  con- 
ducted at  any  number  i of  stress  levels  with  the  constraint,  of  course, 

^ k+1,  allowing  for  estimation  of  the  k+1  polynomial  regression  parameters. 
The  total  number  of  life  tests  conducted  must  be  at  least  8.+1  in  order  that 
at  least  two  tests  be  made  at  one  or  more  stress  levels,  thus  insuring  that 
C is  estimable. 

The  procedure  to  be  described  allows  one  to  obtain  a small-sample  approxi- 
mate prediction  interval,  an  approximate  lower  confidence  bound,  for  the  time 
to  failure  of  the  first  (or  second,  third,  etc.)  failure  in  a finite  lot  of 
specified  size.  To  obtain  the  prediction  interval  one  must  first  calculate 
best  (or  approximately  best)  linear  unbiased  estimates,  best  (or  approximately 
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best)  linear  invariant  estimates,  or  maximum-likelihood  estimates  of  the 
location  and  scale  parameters  of  the  distribution  of  the  logarithm  of  time 
to  failure.  Details  of  the  method  are  described  in  Sections  3 and  4,  pre- 
cision of  the  approximation  Is  described  in  Section  5,  and  an  example  Is 
slven  in  Section  6. 

2.  The  Model 

As  noted  above,  we  assume  that  failure  time  T has  a Welbull  distribu- 
tion with  shape  parameter  C ^ « and  scale  parameter  5 > 0.  Then  X « J.nT 

has  an  extreme-value  distribution  (type-I  distribution  of  the  smallest 
extreme)  given  by 

Fj^(x)  ■ 1 - expj-exp(x-n)/C}  , C > 0 (1.1) 

where  n ■ <in6  is  the  distribution  mode  (a  location  parameter)  and  the  scale 

parameter  is  the  distribution  standard  deviation.  The  parameter  C is 

fixed,  but  n varies  with  level  of  stress  applied  to  induce  failure. 

We  consider  ordered  observables  X,  ,...,X  , r < n,  X..,  > X.  , 

l,n  r,n  i+l,n  l,n 

i = l,...,r,  from  a sample  randomly  selected  from  a distribution  given  by 

(1.1).  That  is,  we  assume  termination  of  the  life  test  at  time  t ••  exr>(x  ). 

r ,n  r ,Ti 

* 

Let  X be  defined  as  the  linear  combination  of  X,  ,...,X  with  smallest 
r,n  l,n  r,n 

mean  squared  error  in  estimating  its  expectation,  of  the  form  r\  + k^, 

-oo  < < < oo.  It  is  thus  the  linear  estimate  most  representative  of  the 
variate  X (logarithm  of  failure  time)  associated  with  the  corresponding  n 
(stress  level). 

Ill  * * ilr 

Mann  [8]  shows  that  X is  equal  to  q - (B/C)^  where  n and  ^ 

2 

are  best  (minimum  variance)  among  linear  unbiased  estimators  of  n and  5,  CL 
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* 2 * -* 
is  the  variance  of  E and  Is  the  covariance  of  n and  C • The 

C^n  ITyll  ITfll 

2 2 2 * 
variance  is  (A-B  /C)C  , where  AC  Is  the  variance  of  n • Note  that  the 
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covariance  of  X and  C is  ~ (B/OC^  = 0.  This  is  a property  that 
r ,n  r ,n 


will  be  Important  in  constructing  our  prediction  interval.  Also  note  that 


for  a complete  sample  from  a normal  distribution  with  mean  U and  standard 
* _ ” 

deviation  a,  X is  equal  to  X = /Z  X. /n.  This  holds  because  the 
n,n  ^ n,n  i 


distribution  is  symmetric  and  hence  the  covariance  of  X^  ^ and  any  unbiased 


linear  estimator  of  a is  zero.  We  comment  here  that  the  results  that  follow 
apply  in  principle  when  X has  a normal  distribution  by  simply  substituting 
y for  n and  O for  C. 

It  is  useful  to  mention  at  this  point  that  tables  of  coefficients,  for 
obtaining  optimal  or  nearly  optimal  linear  estimates  of  extreme-value  or 
normal  parameters  along  with  values  of  A,  B and  C (or  functions  of  these) 
are  given  by  Mann  [6,7],  Mann,  Schafer  and  Slngpurwalla  [14],  Mann  and 
Fertig  [12],  Sarhan  and  Greenberg  [17],  and  Mann,  Witherspoon  and  Fertig  [15] 
as  functions  of  r and  n.  Also,  results  in  Lawless  [4],  Fertig  and  Mann  [3], 
and  Mann  [10]  indicate  that  maximum-likelihood  estimates  can  be  substituted 
for  best  linear  invariant  estimates  in  obtaining  the  approximations  to  be 


described  herein.  Best  linear  invariant  estimators  of  n and  £i  are 


n = n “ BC  and  C “ C /(1+C),  respectively,  so  that 

iTfTi  iTfH  iT^n  iTfH  ir^n 

* ^ ^ 

X = n - (B/C)^  is  equal  to  n - (B/C)s  . This  means  that  the 
r,n  r,n  r,n  r,n  r,n 

* 

best  unbiased  and  best  invariant  estimators  of  E(X  ) ■ n - (B/C)£;  are 

r ,n 


equivalent.  This  is  true  of  no  other  location  parameter  of  the  form  n + <C. 
Note  that  for  the  normal  distribution  best  linear  unbiased  and  best  linear 
Invariant  estimators  of  y + ko  are  not  equivalent  unless  ic  - (-B/C)  = 0. 
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The  model  we  are  assuming  as  relacing  Che  logarlchm  of  failure  Cime 

A 

and  a funcCion  of  scress  level  s is  polynomial.  Since  is  Co  represenc 
Che  logarlchm  of  failure  cime  ac  scress  level  s,  we  have 


E(X*)  - 6q  + 6^1g(s)f  + ...  + 8^jg(s)}’ 


(2.1) 


Using  X^  allows  for  most  efficient  use  of  the  data  in  estimating  or  predicting 

values  of  the  parameter  n + <^  for  any  k or  of  failure  times  at  the  nominal 

level.  See,  for  example,  tabulations  given  by  Mann  [11],  in  which  comparisons 

are  made  of  results  based  on  (2.1)  and  those  based  on 

^ k 

E(X  ) ■ n - yC  ■ X,  6,(g(s)]  with  Y equal  to  Euler's  constant,  approxl- 

® ® s j-0  J 

mately  0.577216. 

The  method  to  be  described  can  be  used  even  though  only  a single  obser- 
vation is  made  at  one  or  more  levels  (as  long  as  the  total  number  of  observa- 
tions i.  ^ k+1) . When  r ■ 1 and  the  sample  size  at  the  ath  stress  level  is  n . 

a’ 

* 

then  X Is  equal  to  X,  with  expectation  r|  - yC  “ ^nn  S,  where  the  sub- 
na  ^ l,n^  ^ a (X  ’ 

xcrlpt  ot  indicates  the  appropriate  value  at  the  ath  level.  The  variance  of 

2 

X,  is  ttC  /6.  Note  that  difficulties  arise  when  any  n is  equal  to  one 

l.n^ 

if  one  attempts  to  use  a model  directly  relating  g(s)  and  p,  as  one  might 
attempt  to  do  for  the  Inverse  power-law  model.  The  difficulties  occur 
because  two  observations  are  needed  at  each  level  to  estimate  n* 


3.  Calculation  of  Approximate  Prediction  Intervals  from  Nominal-Level  Data 

In  [9],  a method  is  described  and  investigated  for  obtaining  from 

Weibull  failure  data  a prediction  interval  for  the  first,  second,  etc., 

failure  in  a finite  lot  of  specified  size.  We  let  Y,  be  the  natural 

1 «in 


A 


logarithm  of  the  time  of  the  tth  failure  In  a lot  of  size  m.  We  assume 

a sample  of  size  n has  been  randomly  selected  from  the  same  population 

of  logarithms  of  failure  times  from  which  the  (sample)  of  size  m has  been 

or  will  be  Independently  selected.  Form  X ■*  n - (B/C)C  > a linear 

r,n  r,n  r,n 

function  of  X,  ,...,X  , all  the  available  ordered  observables  of  the  log 

ifO  r,n 

if  ic 

failure  times,  as  described  In  Section  2.  Tien,  (X  - Y.  )/E(X  - Y.  ) 

X yin  X y IQ 

if 

and  C are  shown  In  [9]  to  have  zero  covariance  and  each  to  be  approximate 
chl-square-over-degrees-of-freedom  variates.  Since  this  is  so,  then  under 
certain  conditions,  relating  the  two  numbers  of  degrees  of  freedom. 


F - (X^  „ - Y )/[C*  „|-B/C  - E(Z  )|] 
y r,n  i,m  r,n'  i,m 


(3.2) 


with  E(Z.  ) ■ E{(Y.  - r|)/Cl.  might  be  expected  to  be  approximately 

i,m  l,m  ’ 

distributed  as  a classical  F variate  with  degrees  of  freedom 


V,  - 2|b/C  + E(Z,  )f^/{A  - B^/C  + Var(Z,  )} 
i l,m  ’ ' l,m  ' 


V,  - 2/C  . 


Values  of  E(Z.  ) and  the  reduced  variance  Var(Z,  ) are  published  by  White 

X yID  X yin 

[19]  for  samples  of  size  m,  m » 2(1)50(5)100.  Also,  if  1 Is  small.  It  is 


relatively  easy  to  calculate  E(Z,  ) from  the  expression  given  bv  Lieblein 

1 ,m 


5 


1 

! 

i 


A similar  expression  In  [5]  Yields  Var(Z.  ).  For  i ■ 1,  E(Z.  ) - -Y-2'nni 

1 im  1 ,m 

2 

and  Var(Z.  ) ■ n /6.  A method  for  calculating  percentiles  of  an  F dlstrl- 
1 ^m 

butlon  with  noninteger  degrees  of  freedom  Is  given  by  Mann  and  Grubbs  [13] 
and  described  In  [14]  and  [9]. 

Combinations  of  values  of  and  for  which  the  F approximation  to 
F^  is  sufficiently  precise  for  use  with  Weibull  or  lognormal  data  is  given 
by  Mann  [9]  and  Fertlg  and  Mann  [3].  Note  that  V2  depends  upon  sample  size 
n and  censoring  number  r while  depends  not  only  upon  r and  n,  but  also 
upon  1 and  lot  size  m. 

To  obtain  an  approximate  lower  prediction  limit  for  ^ at  confidence 
level  1-a  from  (3.2),  one  simply  forms 

where  F^^  2(?(9(l-a)th  percentile  of  the  F distribution  with 

and  degrees  of  freedom.  Because  tends  to  be  small  with  respect  to 
here,  precision  of  results  based  on  the  method  described  in  [13]  are 
greatly  Increased  if  one  approximates  1/F^(v^  ,Vj^) , which  is  equal  to 
Fj^  than  approximating  F,a  directly. 

In  Section  4,  we  describe  how  this  result  can  be  adapted  to  a situation 
in  which  life  tests  are  conducted  at  levels  other  than  nominal. 
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Nominal  Level  Prediction  Intervals 


Let  us  suppose  first,  as  assumed  by  Mann  [8],  that  the  life  tests  are 


to  be  conducted  at  k+1  levels,  where  k specifies  the  degree  of  the  poly- 


nomial relating  E(X^)  and  g(s) . Fixing  the  number  of  levels  at  Icfl 


facilitates  the  exposition  of  the  analysis  by  allowing  for  expression  of 


if  if 

X^,  the  best  linear  unbiased  estimate  of  value  of  E(X  ) at  the  nominal 


level  s^,  in  terms  of  the  Legendre  polynomial. 


(h-hj^)  . . . (h-h^_j^)  . . . (h-hj^) 


a 


(h  -h-  ) . . . (h  -h  , ) (h  -h  . (h  -h,  ) 

a 1 a a-1  a a+1  a Tc 


1,... ,k+l 


where  h » g(s)  and  h^  ■ g(s^).  Subsequently  we  will  treat  the  more  general 


case  in  which  the  number  of  levels  i.  may  be  greater  than  k+1. 
From  Mann  [8],  one  sees  that  under  the  assumption  (2.1), 


k+1 


« VMM  W 

X - y L (h  )X 
o a o a 

ci»l 


(A.l) 


and  that  the  best  (minimum-var lance)  linear  estimator  C of  C Is  given  by 


* ^ 


with  variance  C 


•[i; 


-1 


* * * 
Here  X , C and  C are  the  values  of  X , 
a’  a 


C and  C applying  to  the  sample  size,  censoring  number  and  the  data  at  the 


ath  level.  The  expectation  of  X^  is 


* r 1 

:«„)  - k-  5 E = k - 

L a*i  -* 
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r 


and  the  variance  of  is  given  by 


k+1 

»«(<)  ■ E UaO-o’l'Wc.  = v.r(Z*)5^ 

a-1 


with 


z*  - (X*  - n )/C  . 

o o o 


Icfl 


It  is  Important  to  recall  that  y*  L (h  ) » 1.  This  implies  that 

ct  o,.j. 


k+1 


a»l 


k+1 


X - - T]  L (h  )X  - is  equal  to  Y L (h  ) (X  -Y^  ) and  hence 

o l.m  a o i,m  ^ (W  “ ° 

that  X - Y.  is  a linear  function  of  the  approximate  gamma  variates 
o 1 ,m 

* 

(X  - Y.  ),  a « l,...,k+l.  Therefore,  by  results  of  Box  [17]  and  Patnaik 
a i,m  ■' 

* * 

[16]  (X  - Y,  )/E(X  - Y.  ) is  approximately  a chisquare  over  its  degrees 

o l,m  o l,m 

of  freedom  given  by 


- 2|(B/«^  + E(Z^_j(^/[Var(Z*)  + VarCZ^  j] 


The  estimator  ^ is  also,  clearly,  a linear  function  of  the  s,  each  of 

which  is  approximately  an  independent  chisquare  over  degrees  of  freedom  and 

* * 

5 is  an  unbiased  estimator  of  C-  Thus,  C /?  is  approximately  chisquare 

over  degrees  of  freedom  with  degrees  of  freedom  equal  to 


k+1 

'2  - 2E  . 

a*l 


As  a result,  one  would  expect  that  under  appropriate  conditions  for 


and  v^. 


8 


(4.2) 


Is  approximately  distributed  as  an  F-variate  with  v°  and  degrees  of  free- 
dom. Consequently,  a lower  (l-a)-level  confidence  bound  on  T,  ■ exp(Y.  ) 

i,m  l,m 

is  approximately  given  by 

exp[X*  + + (B/OjC*]  • (4.3) 

If  the  number  i of  testing  levels  is  greater  than  k+1,  then,  the  model 
cannot  be  expressed  as  in  (4.1).  Rather,  we  have, 

E(X*)  - HB  (4,4) 

o — ' 

where  8 is  the  column  vector  (B  ,B, ,---*B,  ) and  H is  a matrix  of  which  the 

^ O 1 iC  ^ 

element  in  the  ctth  row  and  jth  column  is  equal  to  h^  ■ |g(o^)}^. 
a ■ l,...,)l,  j ■ 0,1,..., k.  Then,  the  best  linear  unbiased  estimator  vector 
for  B is 

8*  =»  (hV~^)"^hV“^X*  , (4.5) 

* * * , 

where  X is  the  column  vector  (X^^, . . . .X^^)  and  V is  the  diagonal  matrix  with 

2 * 

elements  consisting  of  (A  ~ B /C  ) , a * 1,...,?..  The  estimator  X is  then 

a a a o 

given  by  the  right  side  of  (4.4)  with  B * (S^.S  . . . ,Bj^)  substituted  for  3. 

The  expectation  -(B/C)^  is  given  by  HW  where  W is  the  vector  defined  by  the 

, , * * 
right  side  of  (4.5),  with  (Bj^/C^^, . . . ,Bj^/Cj^)  substituted  for  (Xj^,  — • 

* ' -1  -12 

The  covariance  matrix  of  the  elements  of  g is  given  by  (HV  {J)  C so  that 

the  variance  of  Z is  h'  (HV  ^H)  ^ , where  h'  ■ (l,g(s  ),...,  ig(s  )i^). 

o ~o  — - ~o  o o I®  of 
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In  the  next  section,  we  give  results  of  a Monte  Carlo  study  undertaken 


to  determine  the  precision  of  the  approximation  (4.2). 

5.  Precision  of  the  Approximation 

A Monte  Carlo  study  revealed  that  constraints  on  the  precision  of  the 
approximation  in  [9]  for  Weibull  failure  data  at  a nominal  level  apply 
roughly  if  we  modify  the  constraint  on  so  that  it  depends  upon  the 
number  I of  stress  levels  at  which  the  life  tests  are  conducted  as  well  as 
upon  sample  sizes  and  censoring  numbers.  Thus,  for  confidence  level  equal 
to  1-ci,  the  approximate  F values  are  in  error  by  less  than  one  percent 
roughly  for 

1-a  - 0.90: 

1-ct  - 0.95: 

1-a  - 0.98: 

a =»  !,...,£ 

All  the  Monte  Carlo  results  herein  have  been  obtained  using  the  approxi- 
mation of  Mann  and  Grubbs  [13]  with  *^2^  obtained  as  [F^(v^,Vj^)l 

as  described  in  Section  4.  When  n » r , a * 1 1,  > 3.88  if  n 4, 

a " 1, ...  ,2.  and  5*  15. 162  if  n^  6,  a ■ 1, ...  ,2. 

Recent  unpublished  results  of  Fertlg  and  Mann  suggest  that  a more 
generally  applicable  approximation  to  the  distributions  of  (3.2)  and  (4.2) 

10 


v°  > 8.882,  0.15v°  + 4.7  < v°  < 1.9v^  - 1.2 
v°  > 8.882,  0.27V2  + 5.6  < v°  < 0.90  + 6.5 
v°  > 15.162,  O.36V2  2.3  < < 0.78V,  + 4.1 


are  given  by  a logarithm  of  an  F distribution  divided  by  an  independent 

chi-squart  variate.  It  is  planned  that  percentiles  of  this  distribution 

will  be  tabulated  for  various  combinations  of  the  three  parameters. 

In  the  15  cases  studl>':d  in  the  Monte  Carlo  analysis,  Y.  was  taken 

1 ,m 

as  and  g(s)  was  taken  to  be  £.n(s)  (applicable  to  the  inverse-power  law 

model).  A Monte  Carlo  sai!g>le  size  of  5000  was  used  and  appeared  adequate 

to  Insure,  in  all  cases,  precision  of  within  about  1 percent  at  1-a  ■ .98, 

It  was  found  that  although  the  approximation  worked  very  well  in 

general  for  k * 1,  i.e. , linear  (inverse  power-law)  models  when  m was 

I 

sufficiently  large  compared  with  ^ n , it  did  not  perform  well  for  k > 1 

a»l 

unless  s^  was  inside  the  range  spanned  by  s^^, . . . ,Sjj^,  For  k > 1 and 

< Sj^  < s^  < . . . < s^,  the  value  of  v°  tends  to  be  too  small  to  satisfy  the 

constraints  given  above.  The  precision  of  the  approximation  increases  as 

the  n *s  increase  for  values  of  p =»  r /n  fixed  and  as  the  p 's  increase 
a *^0  0 0 a 

for  fixed  values  of  the  n 's.  Increase  in  the  value  of  the  lot  size  m,  for 

a 

all  else  fixed,  also  Increases  the  precision  of  the  approximation.  One  would 
also  expect  some  degradation  in  precision  for  increase  in  1 with  fixed  m. 

The  value  of  v?  is  affected  not  only  bv  the  values  of  r and  n , 
a * 1,...,J!.,  and  of  i and  m,  but  also  by  the  relations  of  the  stress  levels 
among  themselves  and  to  the  nominal  level  s^  and  by  the  degree  of  the  poly- 
nomial to  be  fitted. 

Given  in  Table  1 are  results  of  six  of  the  Monte  Carlo  analyses.  In 
each  case,  s^  is  equal  to  5 or  15,  Sj^  is  equal  to  15  or  5,  “ 50  and,  if 

used,  Sj  ■ 100.  As  noted  earlier,  the  function  g(s)  was  taken  to  be  i!.n(s). 
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Table  1.  Approximate  and  Monte  Carlo  (M.C.)  Values  of 

lOOyxM  Percentiles  of  the  Approximate  F-Variate 

F,  : g(s)  - *!.n(s) 
n 
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6.  Example 


Two  random  samples  of  bonded  composite  joints  were  fatigue-to-fallure 

tested  at  accelerated  test  levels.  A Weibull  distribution  is  known  to  be 

appropriate  as  a failure  model.  A sample  of  size  7 was  subjected  to  test 

* 

at  stress  level  s^^  ■ 20  until  four  failures  were  observed,  and  y ■ 5.51 
* 

and  y • .064  were  calculated  from  the  observed  number  of  cycles  to 
failure.  A second  sample  of  size  5 tested  until  all  five  specimens  failed 

if  if 

at  s,  ■ 50,  yielded  n,.  c ■ 4.31  and  Cc  c * .058.  The  Inverse  power-law 
model  was  used  to  predict  the  time  of  the  first  failure  in  a lot  of  size 
100  at  the  nominal  level  s^  - 10. 

From  tabulations  in  [14],  we  find  for  r^^  ■ 4,  n^  ■ 7 that  A ■ 0.37265, 

B ■ 0.15696  and  C ■ 0.28016  so  that  X,  _ * 5.47.  Similarly,  from  A ■ 0.23140, 

it 


B - -0.033991  and  C • 0.16665  for  r^  - 

Since  hg  ■ UnlO  - 2.302,  h^^  ■ i?.n20 

I /K  \ 2.302  - 3.912  , ,CQ  . T /■u 
^l^^O^  * 2.996  - 3.912  " 1*758  and  L2(h 

X*  - 1.758(5.47)  - 0.758(4.32)  - 6.34. 


“ 5 we  have  X^  ^ * 4.32. 

2.996  and  h2  * )!.n50  ■ 3.912,  then 


N . 2.302  - 2.996 
0-’  ” 3.912  - 2.996 
* 


-0.758.  Thus, 


Also , ^ 


a,  ,/0. 28016  + ,/0. 16665)/ 

4,7  5,5 


1 


(1/0.28016  + 1/0.16664)  - .0602  with  variance  0.10449^“. 

* 


The  expectation  of  (Xq  - Is  then  -[1.758(0.56025)  - 

0. 758(-0. 20398)  -0.577216  - 4.605]  - 4.043  and  the  variance  of  (Xq  - 
is  1.758^(0.37265  - 0.08794)  + 0.758^(0.23140  - 0.00693)  + 1.6449  - 2.653^. 
Thus,  the  degrees  of  freedom  for  the  approximate  F-varlate  F^^  are 
u°  ■ 2(4. 043^) /2 .653  ■ 12.33  and  v°  ■ 2/0.10449  ■ 19.14.  Using  the  method  of 
Mann  and  Grubbs  [13]  we  approximate  first  the  tenth  percentile  of  a Beta 
distribution  with  parameters  19.14/2  - 9.57  and  12.32/2  » 6.16.  This  is 
accomplished  by  letting  m equal  iln(9.57  + 6.16  ■ 0.5)  - in(9.57  - 0.5)  » 0.5183 
and  V equal  (9.57  - 0.5)"^  - (9.57  -*•  6.16  - 0.5)“^  - 0.4493  in 
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0.4993 


I 


Then  F ^^(19. 14, 12. 32)  (12.32/19.14)0.4493/(1  - 0.4493)  - 0.525  and  j 

F gQ(12. 32,19. 14)  - (0.525)“^  - 1.903.  ! 

To  calculate  a prediction  interval  by  use  of  (4.3)  for  the  first 
failure  in  a lot  of  size  100,  we  need  the  values  of 

<®/c>0  ■ -«’'o  - ’'l,100>'5  • • 

i 

1 

* * 

Z “ 0.0602  and  Xq  “ 6.34.  Then,  from  (4.3),  a 90  percent  prediction  inter- 
val for  the  first  failure  in  a lot  of  100  at  the  nominal  stress  level, 

Sq  ■ 10,  is  given  by 


exp[6.34  + 1.903(-4. 043)0. 0602]  - exp(5.88)  - 357 


numbers  of  cycles  to  failure. 

Since  and  easily  satisfy  the  conditions  specified  in  Section  5 
for  precision  of  the  approximation  to  be  within  1 percent,  the  approximate 
prediction  interval  for  should  be  correct  to  within  a unit  In  the 

second  decimal  place  and  the  approximate  prediction  interval  for  exp(Y 
should  be  within  1 percent  of  the  exact  value. 
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